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ABSTRACT 

We study microlensing of the faint images that form close to the centers of strong 
gravitational lens galaxies. These central images, which have finally begun to yield to 
observations, naturally appear in dense stellar fields and may be particularly sensi- 
tive to fine granularity in the mass distribution. The microlensing magnification maps 
for overfocussed (i.e., demagnified) images differ strikingly from those for magnified 
images. In particular, the familiar "fold" and "cusp" features of maps for magnified 
images are only present for certain values of the fraction / of the surface mass density 
contained in stars. For central images, the dispersion in microlensing magnifications 
is generally larger than for normal (minimum and saddle) images, especially when 
the source is comparable to or larger than the stellar Einstein radius. The disper- 
sion depends in a complicated way on /; this behaviour may hold the key to us- 
ing microlensing as a probe of the relative densities of stars and dark matter in the 
cores of distant galaxies. Quantitatively, we predict that the central image C in PMN 
J1632— 0033 has a magnification dispersion of 0.6 magnitudes for Rs,-c/Re ^ 1, or 0.3 
mag for Rsrc/RE = 10. For comparison, the dispersions are 0.5-0.6 mag for image B 
and 0.05-0.1 mag for image A, if Rgrc/RE ^ 1; and just 0.1 mag for B and 0.008 mag 
for A if i?src/-RE = 10. (The dispersions can be extrapolated to larger sources sizes 
as fj cx R^^l-) Thus, central images are more susceptible than other lensed images to 
microlensing and hence good probes for measuring source sizes. 

Key Vifords: galaxies — stellar populations — gravitational lensing 



1 INTRODUCTION 



Two aspects of strong gravitational lensing that date 
back 25 years have attracted much recent interest. 
Microlensing, or variations in the (optical and/or 
x-ray) fluxes of lensed quasar images that are in- 
duced by individu al stars in the lens galaxy, was 
first discussed by IChang fc Refsdall (|l979l '). The phe- 
nomenon has now been observed in several lens systems 
JWozniak et al. 2000: Sch cchter et al. 2003^ Richards et al 
20041: lGavnullina_et all boosl : iKochanek et al.l l2006l : 

be used to learn about 



Morgan et all 120061 '). It 



the st e llar components llGott 198ll: Wambsganss et al.l 
1990b; Pelt ot al. 'l998'; Schmidt fc Wambsganss! Il998l : 
Wyithc ct al. 2000; Kochanck 20 04) and the relative den- 
sities of sta rs and dark matter ijSchechter fc Wambsgans3 
|2002| . |2004| '1 in lens galaxies, and to probe structure in 
the sou rce quasars at stu n ning micro- arcsecond reso- 
lution l|Grieger et all 1 19881 . Il99ll : lAgol fc Kroli^ 1 19991 : 



Mineshige fc Yoneharal Il999l: iFluke fc Websterl Il999l : 
Richards et al.ll2004 iGoicoechea et al.ll2004h . 



Centr al or "odd" lensed im ages were o rigina lly pre- 
dicted bv lDver fc Roedeij l|l980l) and iBurkd (|l98ll ). If the 
inner surface mass density of a lens galaxy is shallower than 
E oc and there is no point mass at the center, then lens- 
ing should always produce an odd number of images. Nearly 
all observed lenses, however, exhibit two or four images. 
The prediction and observations can be reconciled by noting 
that one of the expected images should be very close to the 
center of the lens galaxy and demagnified by the high cen- 
tral surface density t here fe.g..[Walling:ton fc NarayarJIiggSl : 
iNorburv et al. I I2OOOI : iRusin fc Mall200ll : lKeetonll2003l ). That 
makes central images difficult to det ect — and indeed the 
first ones have only just been found. W inn et al. I ([2003, 200i) 
detected a central image at radio wavelengths i n the asym- 
metric 2-image lens PMN J1632— 0033, while llnada et al.l 
l|2005l ) detected a central image at optical wavelengths in 
the unusual 4-image cluster lens SDSS J1004-I-4112. 
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Image 




7 




A 


0.370 


0.422 


4.570 


B 


3.309 


2.859 


0.352 


C 


16.517 


13.796 


0.020 



/ cfl cff\ 

[i^ ,7 



(3) 



Table 1. Total convergence k, shear 7, and magnification n for 
each of the three images in PMN J1632— 00 33, from spherical 
powe r law plus external shear lens models by IWinn et al.l (|2003| . 
l2004h . "C" indicates the central, highly demagnified image. 



Our goal in this paper is to combine the two phenomena 
and study microlensing of central images. The issue is timely 
since central images are now being observed, and compelling 
since the images naturally form in regions where the density 
of stars is high and microlensing seems inevitable. 

We customize our calculations to the lens PMN 
J1632— 0033, because it is both a known central image sys- 
tem and also a prototype for the sorts of systems that 
are expected to y ield the most centra l images in the fu- 
ture (Kccton 2003!: lBowman et aL|[200i ). IWinn et al.l l|2003l . 
12004 ') recently studied a wide range of mass models for PMN 
J1632— 0033, which can be used to estimate the total con- 
vergence K and shear 7 at the positions of the images. The 
lens data, including the central image, are consistent with 
a simple power law surface mass density E cx with 
a = 0.91 ± 0.02, plus a small external tidal shear. Table [1] 
gives the values of « and 7 for such a model, which we adopt 
as inputs for our microlensing simulations. 



2 METHODS 

We use standard methods for microlensing calculations 
(|Wambsgansslll999l ). Specifically, we pick a patch around an 
image that is large compared with a stellar Einstein radius 
(the scale for microlensing), but small compared with the 
global scale of the lens so that the mean densities of stars 
and dark matter are essentially constant across the patch. 
Attaining this balance is not difficult because the global scale 
is set by the l'.'5 image separation while the stellar Einstein 
radius is i?E ^ 2 (M/Mq)^^^ x 10~® arcsec. In practice, we 
typically consider a rectangular patch in the image plane 
chosen so that it maps into a square patch in the source 
plane that is WORe on a side (see below). Without loss of 
generality, we can choose a coordinate system aligned with 
the direction of the local shear. We assume that a fraction 

/ = ^stars/'^tot (1) 

of the total surface mass density (or convergence) is con- 
tributed by stars. We use that to determine the number 
density of stars in the image plane patch, 

/ Ktot 



(2) 



and then distribute the stars randomly. The total surface 
density Ktot is kept fixed (to the values in Table [l]) by 
including a sm ooth, continuous matter component with 
Kc = {1— f)Ktnt ■ [Paczvnskil (|l986l ) has shown that including 
a smooth matter component is equivalent to rescaling the 
convergence and shear to 



with / — 100% in stars. For the case of a finite source, this 
also leads to a rescaling of the source size, 

^src 



(4) 



In this pilot study, we assume that all stars have the 
same mass, and we always work in units of the stellar Ein- 
stein radius _Re. Previous analyses have shown that the mi- 
crolensing magnification distribution is at best weakly de- 
pendent on the distribution of stellar masses, at least for 
a point sourceQ That result must break down for a finite 
source, which would be insensitive to stars below some mass 
threshold (roughly corresponding to _Re ^ Rsrc)- The full 
problem — microlensing of a finite source by stars with 
unequal masses — is c ertainly interesting, b ut beyond the 



scope of this paper fcf. ICongd on et al 



We use ray-shooting softw are by 



20061) 



Wambsgansj l|l990D . 



IWambsganss et all (|l99()lj ). and lWambsganssI l| 19991 ) to per 

form the microlensing simulations. Briefly, the software 
"shoots" a large number (> 10*) of light rays from the ob- 
server through the lens/image plane into the source plane, 
and coUects them in small pixels. The number of rays in 
a source pixel is proportional to the lensing magnific ation 
at that position. As described by Wambsganss (1999), the 
rectangular shooting patch in the image plane is larger than 
the region that would map into the source plane box in the 
absence of microlensing. The reason is that the defiection 
diverges as 1/r near a star, so stars outside the nominal 
shoo ting region can s end light rays into the source plane 
box (|Katz et al.lll98d ). The shooting region is expanded to 
ensure that a significant fraction of this "diffuse flux" is col- 
lected. 

The ray-shooting software produces a magniflcation 
map in the source plane. We can obtain the magniflcation 
probability distribution by making a histogram of the pixel 
magniflcations. To consider a source with a flnite extent, we 
convolve the magniflcation map with the surface brightness 
distribution of the source before making the histogram. For 
simplicity, we use a Gaussian sour ce with half-light radius 
Rbtc, since [M"ortonson et al.l ((200^ argue that the detailed 
structure of the source hardly affects magniflcation distri- 
bution s provided that th e half-light radius is the same (also 
see Congdon et al.l |2006|). 

The physical input parameters for the simulations are 
the total convergence and shear, which we take from Table 
[T] and the fraction / of mass in stars, which we vary (see §3). 
The technical input parameters are the size of the magnifi- 
cation map and the number of pixels. We seek to consider 
sources both smaller and larger than _Re, and in particular 
find the range 0.1 ^ Rsrc /Re ^ 10 to be both interesting 
and tractable. In order to have many independent source po- 
sitions in a given magnification map, we use maps that are 
WORe on a side. In order to have enough pixels to handle 
small sources, we use maps with 1024 x 1024 pixels. 



^ The magnification distribution w as long thought to be indepen- 
dent o f the microlens mass f u nction | Schneid eilll987t|Wambsganssl 
1992"; 'Lewis & IrwinI Il995l : IWvithe fc TurneJ I2OOII) . However. 
Schochtcr ot al. ( 2004i) recently presented results from numerical 
experiments suggesting that there is a weak dependence. 
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To obtain fair sampling, we create 50 independent mag- 
nification maps for each image, and compute the final mag- 
nification distributions from the combination of all of them. 
A fringe benefit of this approach is that we can consider 
the 50 magnification maps to be independent realizations of 
the microlensing calculation, and use bootstrap or jackknife 
resampling (e.g.. lEfronlllQsd ) to estimate the statistical un- 
certainties in our analysis. 



3 RESULTS 

3.1 Magnification maps 

Figure [T] shows magnification maps for / — 100% of the 
surface density in stars. To display the range of features 
that appear on different length scales, we show (1QQ_Re)^, 
(30-Re)^, and (10-Re)^ maps here. The maps for image C in 
particular show that there is structure on scales of tens of 
Einstein radii, so even a (100-Re)^ map does not contain a 
fully representative sample of magnifications. 

The qualitative differences in the three magnification 
maps shown in Figure [1] are striking. For image A, which 
forms at a local minimum in the time delay surface and 
has relatively low values of the convergence and shear lead- 
ing to a modest ampli fication, the map shows the famil- 
iar caustic network (e.g. . I WambsgansJiggol : ISchneider et al.l 
1 19921 : IWambsgansd Il999i r The caustics are preferentially 
stretched along the horizontal axis because we have cho- 
sen coordinates aligned with the local shear. For image B, 
which forms at a saddle point in the time delay surface and 
is modestly demagnified, careful inspection reveals that the 
map shows many of the 3-pointed cusps that are charac- 
teristic of microlensing when the parity is negative (e.g., 
ISchneider et al.ll 19921 : [Petters et aLlboOlj ). These resuhs are 
familiar. More surprising is the map for image C, which 
forms at a local maximum in the time delay surface and 
is highly demagnified. This map looks very different from 
the familiar, traditional caustics. It exhibits round, con- 
centrated blobs of high magnification, surrounded by large 
regions of low magnification. These features are known to 
occur for highly o verfocussed (i.e., demagnified) microlens- 
ing scenarios (e.g., 'Wambsganss' 'l99rt JSchneider et al.lll992l : 
iPetters fc Wittlll996 : Petters et al..,2Q()ih . but thev have not 
been explored in great detailQ 

As noted above, there appears to be strong inhomogene- 
ity even on scales of tens of Einstein radii. To quantify this, 
we define the auto-co rrelation function along l inear slices 
through the maps (see I Wambsganss et al. I ll990al ): 



{^l{x)^l{x + Ax)) - (^Lix)f 
{ti{xY) - {t,{x)y 



(5) 



where fi{x) is the magnification at position x and the aver- 
ages are over all positions along the slice. For each map, we 
calculate £,{Ax) for each row and average over rows. We then 
further average over 50 maps (i.e., 50 random stellar configu- 
rations). This gives ^(Ax) parallel to the shear axis; we then 



2 IPetters fc WittI l|l996t ) and lPetters et al.l l|200ll . see Chapter 15) 
have shown that, in the case of an isolated star with a constant 
shear and a large, constant convergence Kc, the caustics consist 
entirely of elliptical fold curves and do not have any cusps. 



repeat the procedure averaging over columns to find ^{Ax) 
perpendicular to the shear. The auto-correlation function for 
each image, both parallel and perpendicular to the shear, is 
shown in Figured The distance Axq at which ^(Ax) falls 
to zero is a measure of the correlation length. 

Parallel to the shear, Axq ~ 6.1-Re while Axq « 
10.67?E, implying that the image B maps have structures 
spanning ~10-15 -Re and the image C maps have structures 
~20-25 -Re. It is difficult to interpret parallel to the 
shear since the caustics are preferentially stretched along 
the shear axis, leading to substantial power in even at 
large Ax. Perpendicular to the shear, falls rapidly to 
zero {Axq « 2.3-Re along this direction) and in fact Figure 
[1] shows little in the way of large scale structure in the ver- 
tical direction. For images C and B, we see that > for 
all distances Ax. This supports the qualitative conclusion 
that there is more structure on large scales for the image C 
compared to image B maps. 

So far we have assumed that / = 100% of the sur- 
face mass density for each image is in stars, which is prob- 
ably valid for image C but unlikely for image A. In Fig- 
ures [SHS] we examine the magnification maps for differ- 
ent values of /. In Figure |3] when / is low, the map ex- 
hibits caustics of individual stars, and as / increases the 
number of caustics increases a nd they begin to merge (cf. 
ISchechter fc Wambsganss| |2002l ). 

Figures|3]and[S]show drastically different behaviour. For 
image C (Figure (5)1, the caustics morph from blobs to cusps 
as / decreases; by / = 0.2 the magnification map consists en- 
tirely of fold and cusp caustics. Decreasin g / further causes 
holes to appear in the magnification map (Chan g fc Refsdal 
119841 . For image B, we find the curious result that cusp 
caustics are present for / = 1.0 (Figure [T]), but they are 
absent for 0.9 > / > 0.5, and then they reappear at lower / 
val ues. This behaviour is surprising because it was not seen 
by ISchechter fc Wam bsgansg (|2002l ) in their magnification 
maps for a highly magnified saddle image. Since our im- 
age B is a demagnified saddle image, we conjecture that the 
overall magnification, in addition to the parity, is important 
for determining whether the caustics are blobs or cusps and 
folds. The dependence on the stellar mass fraction appears 
to be complicated. 

3.2 Magnification distributions 

The magnification distributions for different source sizes and 
different images are shown in Figure[6] (assuming / — 100%). 
The magnification is expressed as the change in magnitude 
due to the presence of microlensing. Am = — 2.5 log(/i//io)- 
In general, the distributions appear to be smooth and fairly 
symmetric about Am — 0. The exception is the case of 
image A with a small source {Rsic/Re ^ 0.1), whose magni- 
fication distribution shows three distinct peaks. The peaks 
can be understood in terms of the caustic network (e.g., 
Ranch et al.l 1 19921 : IWambsganss et ah 1 Il992l : iGranot et al.l 
20031 ): from high to low Am — i.e., low to high fj, — , the 
first peak corresponds to sources outside all caustics, the 
second peak to sources inside a single caustic (so there is 
one additional pair of microimages), and the third peak to 
sources inside a pair of overlapping caustics (so there are 
two extra microimage pairs). In fact, there are higher-order 
peaks corresponding to source positions where more caustics 
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Figure 1. Sample magnification maps for images A (left), B (center), and C (right). The top row shows 100-Re on a side, the middle 
row shows 30-Re on a side, and the bottom row shows 10-Re on ^ side. The zoom in region in each case is centered on the (100-Re)^ 
map. The fraction of surface mass density in stars is / = 100% for all cases, and the colorbar at the bottom indicates the change in 
magnitude, with respect to the mean magnification, at each source position (yellow regions indicate magnification, blue regions indicate 
demagnification) . 



overlap, but those peaks blur together and do not appear 
as clear features in the final magnification distribution (see 
iGranot et al.|[2003l . especially their Fig. 4). The multi-peak 
structure disappears as the source size increases, because a 
large source {Rare /Re ^ 1) generally extends over one or 
more caustics, so the regions of different image multiplicity 
in the source plane are smeared out and disappear as dis- 
tinct maxima in the magnification distribution. The multi- 
peak structure is not apparent in the histograms for images 
B and C, even for sources as small as _Rsrc ~ 0.1-Re, because 



the caustics are so small that only a small fraction of source 
positions lead to additional microimage pairs. 

An important quafitative feature of the magnification 
distributions is that the width decreases as the source size 
increases. The width can be considered to represent either 
the fikelihood that the magnification at any given time is 
different from what a smooth model would predict, or the 
RMS amplitude of variations in time. Thus, we may say that 
the larger the source, the less it is affected by microlensing — 
a well known result which makes intuitive sense. To quantify 
this effect, we compute the dispersion in magnitudes, namely 
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Figure 2. One-dimensional auto-correlation functions for magnification maps with / = fOO%. The left panel shows ^(Aa::) aligned with 
the shear axis, while the right panel shows ^{Ax) perpendicular to the shear axis. In the left panel, there is significant power in 
because the caustics tend to be stretched along the shear axis. The fact that > implies that there is structure in the image C 
maps on larger length scales compared to the image B maps. 



a= [((-2.51og^i)2)-(-2.51ogM) 



211/2 



(6) 



where the average is over source positions. In previ- 
ous work the magnification dispersion has been stud- 
ied analyticall y (for the case of 7 = 0, «e S, and 

J^src > 5J^e; iRefsdal fc Stabelll Il99ll: ISeitz et all 
Neindorj 20031') , nuraerically (jPeguchi fc Watsonl 



1994 



1987 



Refsdal fc Stabelllll99ll: [Seitz et al.lll994l : iRefsdal fc Stabelll 
19971 : IWvithe fc Tu rner 2002^), an d with ray-shooti ng meth- 
ods (ISchechte r fc Wa mbsganssi I2OO2I: [Wyithe fc Turneij 



20021: iLewi s fc Ibataj 



jsganssi 
I2OO4I : 



Schechter et al.l 



2004 



Mortonson et al.ii2005 '). 

Figure [7] shows that the dispersion falls monotoni- 
cally as the source size increases (as seen in the previ- 
ous work). It remains relatively constant until the source 
size becomes comparable to the stellar Einstein radius, and 
then drops rapidly. Curiously, all three images have roughly 
the same magnification dispersion when the source is small 
{Rsrc/RK < 0.3). For Rsrc/RE > 0.3, the ordering of the im- 
ages is independent of source size: image C is most affected 
by microlensing, followed by image B, and finally image A. 
For large sources analytic arguments suggest that the disper- 
sion should scale as cr oc 

R^rl (Refsdal & StabcU 1991, 1997,). 
We certainly see this for image A. It appears that images 
B and C are approaching the asymptotic scaling but have 
not quite reached it at Rstc/Re ~ 10. This is not surprising 
because the predicted scaling applies when the number of 
stars in front of the image is large. Since images B and C 
are demagnified, the source needs to be larger (compared 
with the magnified image A) in order for the image to be 
large enough to intercept a significant number of stars. 

The next step is to consider how the magnification 



1.25 



0.5 



E 0.25 



0.125 



0.05 















'■•..A 







0.1 



10 



Figure 7. Magnification dispersion as a function of source size, 
for the three images, assuming / = 100% of the mass in stars. 
The line segment below the curve for image A s hows the expected 
asymptotic scaling for large sources, a oc R~l llRefsdal fc Stabelll 
Il99ll.1l997t ). 



dispersion changes as we decrease /, as shown in Figure 
[S] For images A and B, the dispersion falls monotoni- 
cally as ,f decreases. For image B, this is in contrast to 
ISchechter fc Wambsgans3 (|2002l ) who found that decreas- 
ing / can increase the magnification dispersion for saddle 
images. We believe the difference is explained by the fact 





Figure 3. Magnification maps for image A with size (SORe)^ and varying / (cf. Figure [T] left column, for / = 1.0) The visual contrast 
in each map conveys the magnification dispersion. The dispersion increases with increasing /: rapidly for small /, more slowly above 
/ > 0.3 (compare to Figure |8] with -Rsrc = 0.1). 



that ISchechter fc Wambsgans^ considered a highly magni- 
fied saddle image {fi = 9.5), such as would be seen in a 
bright minimum/saddle pair straddling a critical curve in a 
4-image lens; whereas our saddle image B has a modest de- 
magnification {fi = 0.35), as is common for the image nearer 
the lens galaxy in an asymmetric 2-image lens. 

For the central image C, we find that decreasing the 
fractional surface density in stars can increase the magni- 
fication dispersion, so long as iisrc ^ Re- Specifically, for 
a large source Rare /Re ~ 10 (right panel of Figure O, the 
magnification dispersion for image C decreases monotoni- 
cally as / decreases, like it does for images A and B. When 



Rstc/Re = 1, however, the dispersion increases as / de- 
creases all the way down to / f» 0.04, and only then turns 
over. (The turnover is reassuring, because for / — > there is 
no microlensing and the magnification dispersion must van- 
ish.) When Rsic/Re = 0.1, the turnover occurs at / < 0.01. 
While the behaviour of image C at low / may have little 
practical importance, it is still useful for understanding how 
the magnification dispersion depends on the relative densi- 
ties of stars and dark matter. 

We note that, in the limit of an infinitesimal source, the 
dispersion is visually apparent in the contrast of the magni- 
fication maps shown in Figures |3]|4l and [5] The dispersion in 
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Figure 4. The same as Figure |3] but for image B (cf. Figure [T] middle column, for / = 1.0). The relative number of cusps and folds 
versus blob caustics has a complicated dependence on /. The visual contrast in each map conveys the magnification dispersion. The 
dispersion increases with increasing /: rapidly for small /, more slowly above / > 0.3 (compare to Figure[8]with Rbtc = 0.1). 



the maps of images A and B visibly increases with /, while 
the dispersion of the image C maps decreases with increas- 
ing /. Further, the contrast changes sharply at low / for all 
three images and then remains roughly constant for / > 0.4. 
This is consistent with the slopes of the a versus / curves 
shown in Figure ID 

To make a final interpretation of Figure |S] it is useful to 
estimate an upper bound on the fractional surface density in 
stars at the positio n of each image. For PMN J1632— 0033, 
|w inn et al.] l|2003l ) show that the lens data are consistent 
with a total surface mass density of the form Ktot oc _R~" 
with a — 0.91 ± 0.02. They note that the lens galaxy is too 



faint to be characterized in detail with existing HST images, 
but appears to be an early-type galaxy with effective radius 
of i?e ~ 0'.'2. To estimate the maximum possible stellar mass 
density as a function of radius, we model Katars using a de 
Vaucouleurs R^^''' law profile, and set the stellar mass-to- 
light ratio to the largest allowed value such that Kstars never 
exceeds Ktot- Using the observed distances to the center of 
the lens (Ra = l'.'38, Rb = 0'.'09, and Rc = O'.'Ol; see 
I Winn et al.ll2003l ). we obtain the /max values listed in Table 

m 

This analysis confirms our guess that the matter at im- 
age C could be essentially all stars. By contrast, the galaxy 
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Figure 5. The same as Figure[3] but for image C (cf. Figure[T] right column, for / = 1.0). The relative number of cusps and folds versus 
blob caustics has a complicated dependence on /. The visual contrast in each map conveys the magnification dispersion. In contrast 
to Figures |3] and |4] the dispersion decreases with /, rapidly at first and then more slowly above / > 0.4 (compare to Figure |8] with 



0.1). 



appears to be concentrated enough that at image B the den- 
sity is no more than half stars, and at the distance corre- 
sponding to image A the density is no more than 1% stars. 
When we combine the upper hmits on / with the cr versus 
/ curves in Figure |8l we obtain the upper bounds on the 
magnification dispersion hsted in Table [2] 



The bottom line is that microlensing is basically negli- 
gible for image A, but reasonably important for both images 
B and C. For large source sizes, the central image C is no- 
tably more affected by microlensing than the saddle image 



B, hence microlensing of such demagnified images can help 
reveal the source size. 



4 DISCUSSION 

Central lensed images, which are highly demagnified and 
naturally appear in places where the density of stars is high, 
are more susceptible to gravitational microlensing than the 
more familiar images that form at larger distances from lens 
galaxies. Therefore, it is useful to understand microlensing 
of central images in detail. 
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Figure 6. Magnification distributions for the three images (left to right), for three different source sizes (top to bottom), assuming 
/ = 100% of the mass in stars. The source is modeled as a Gaussian with half-light radius -Rare (quoted in units of -Re)- The distributions 
for each image are derived from 50 independent realizations of (100-Re)^ magnification maps. Each probability distribution P{logfi) 
is normalized to unit area, but the vertical axis scale varies in the bottom panels. Any apparent jaggedness is numerical noise. The 
double peak in the panel at the top left is real and reflects the first additional microimage pair (for details see text). The histograms are 
centered on Am = — 2.5 log(/i//io) = Oi where fio is the magnification in the absence of microlensing (from Table [T]) The microlensing 
distributions yield a mean magnification {/i) that agrees with /iq to within numerical uncertainties. 



The microlensing magnification maps for demagnified 
images difi'er qualitatively from those for familiar magnified 
images in striking ways. For high fractional surface density 
in stars, the central image maps do not show any of the 
classic fold and cusp caustic structures, but rather have con- 
centrated blobs of high magnification amid large regions of 
demagnification. As / is decreased the maps begin to de- 
velop more and more cusps and folds. For the demagnified 
saddle, low values of / yield cusps and folds, intermediate 
values produce blobs, and at very large values (/ « 100%) 
the cusps and folds reappear. Thus the relative number of 
cusps and folds versus blo bs has a co mplicated dependence 
on Ktot, 7, and / (see also lPetters et al..,2001. ). 



The central image C maps also show notable inhomo- 
geneity on scales of tens of Einstein radii. These large scale 
structures are evident in the 1— D auto-correlation func- 
tion ^(Ax). For the central image C, falls to zero by 
Aa;g ~ 10.6-Re along the direction parallel to the shear, im- 
plying structure sizes ~ 20-25-Re. In contrast Axq ~ 6.1 
suggesting smaller structures with sizes ~ 10-15-Re. The 
larger scale structures in the image C maps causes microlens- 
ing fluctuations to be larger for demagnified central images 
than for other images, even when the source is fairly large 

{Rsrc/RE > 1). 

Our most intriguing qualitative result concerns the 
sensitivity of central image microlensing to the relative 
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Figure 8. Magnification dispersion as a function of the fractional surface density in stars, / = Kstars/Ktot, for three different source 
sizes (quoted in units of R^). The three Unes show the three images A (bottom), B (middle), and C (top) of PMN J1632— 0033. Note the 
widely diffe rent scales for a. The errorbars indicate the statistical uncertainties estimated with bootstrap and jackknife resampling (e.g., 
lEfronlllQsj) . 



Image 


/max 


(T (in 

Rsrc =0.1 


mag) 
1 


10 


A 


0.013 ±0.001 


0.13 


0.05 


0.008 


B 


0.511 ±0.023 


0.58 


0.48 


0.10 


C 


0.997 ± 0.003 


0.63 


0.58 


0.28 



Table 2. Column 2 gives our estimate of the upper limit on / = 
Kstars/^tot (see text). The quoted uncertainties are derived from 
the uncertainty in the powe r law index of th e lens galaxy mass 
distribution {a = 0.91±0.02: lwinn et al.l2003l) . Columns 3-5 give 
upper limits on the magnification dispersion cr (in magnitudes), 
for three source sizes (quoted in units of Re)- These upper limits 
were obtained by combining /max with the cr vs. / curves in Figure 
[8] The uncertainties in a due to the uncertainties in /max are 
^ 0.13 mag. The dispersion can be extrapolated to larger sou rce 
sizes with the scaling ct oc R~l llRefsdal fc Stabell|[l99ll . Il997l ). 



densities of stars and dark matter. When the source is 
large {Rbic/Re 3), the magnification dispersion decreases 
monotonically witii tiie fraction / of density in stars. How- 
ever, when tiie source is small {Rstc/Re ^ 3), the mag- 
nification dispersion rises as / is decreased from unity, 
reaches a peak at some finite value of /, and then falls to 
zero as / — > (as it should in the absence of microlens- 
ing). This dependence is simila r to the behaviour seen by 
ISchechter fc Wambsgan"s3 l|2002l ) for a highly magnified sad- 
dle image. We do not see such behaviour for a demagnified 
saddle, but we do see it for a demagnified central image. 
While this qualitative result has little practical importance 
(we expect f ~ 1 for central images), it may be valuable 
as insight into how microlensing depends on the relative 
amounts of stars and dark matter. If we can develop a full 
understanding of that problem, we may be able to turn mi- 
crolensing into the best tool available for probing local densi- 
ties, rather t han integrated masses, of dark matter i n distant 
galaxies (see ISchechter fc Wambsgans3l2002l . |2004| ) . 

One practical goal of microlensing studies is to probe 
the structure of the optical continuum emission regions 
of quasars at very high effective spatial resolution. The 



light curves produced by relative motion of the lens 
galaxy and source quasar can be used to map out 
the (1-dimension al) structure of the quasar on micro- 
arcsecond scales llGrieeer et aLlll988l. Il99ll: lAgol fc KrolikI 



19991: iMineshige fc YoneharaHl999l : iFluke fc Webstedll999l : 
Goicoechea et al. 1 12004 ). In terms of raw variability ampli- 
tude, it would seem that central images are the best targets 
for microlensing observations. Of course, there is the prob- 
lem of detecting central images at optical wavelengths in the 
first place. A faint central image may be swamped by light 
from the lens galaxy. It may also be dimmed by extinction 
or scattering in the interstellar medium of the lens galaxy, 
although those effects might not be too much of a concern 
in the vast majority of lens galaxies that are ellipticals. 

It is interesting to consider the possibility that central 
image microlensing may allow us to probe the size of the 
radio emission region of the source quasar. For large sources 
the dispersion scales as a oc R^^l, so for radio sources with 
Rstc/Re <; 100 we expect microlensing fluctuations to be no 
more than a few percent|f| For the case of PMN J1632-0033, 
the intensity of the central image is -^7 50 and 500 ^Jy at 8 
and 22 GHz, respectively l|Winn et al.li2004 ), so microlens- 
ing fluctuations may be at the level of tens of /iJy. Presently, 
integration times for observing these fluctuations are pro- 
hibitively long for a monitoring campaign (~5 or 50 hrs at 
8 or 22 GHz to resolve 5% fluctuations However, next 
generation telescopes such as the EVLA or EMERLIN will 
make these observations feasible. 

Our results allow us to consider how microlensing may 
affect various other applications that involve central lensed 
images. In PMN J1632— 0033, the position and brightness 
of the central image lead to strong constraints on the den- 
sity profile of the lens galaxy, and to upper limits on the 



However, iKoopmans fc de Bruvnl l|2000h report that there is 
evidence for microlensing fluctuations in at least one radio lens, 
B1600+434. Additionally, scintillation from the Milky Way may 
still have an appreciable effect at these source sizes for low galactic 
latitudes (e.g. Koopm ans et al. 2003 ). 
* c.f., |http://www.vla.nr ao.edu/astro7] 
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mass of any supermassive black hole that may resid e at the 
center of the lens galaxy l|Winn et all |2003| . 120041 '). Those 
constraints are based on radio data where, as noted above, 
the microlensing fluctuations are on the order of a few per- 
cent. However, even if they were as large as tens of percent, 
that would still be negligible compared with the orders of 
magnitude over which the macro magnifications of central 
images can vary d ue to modest changes in the smooth lens 
model (e.g.. lKeeton 2003 1. In other words, microlensing does 
not appear to be a significant concern for constraints on lens 
models drawn from radio observations of central images. 

Perhaps even more interesting is the possibility that any 
central macro image produced by a lens galaxy containing 
a supermassive black hole should be accompanied by a sec- 
ond, fainter central image l|Mao et al.ll200lt iBowman et al.l 
[2OOJ), and that the central image pair c an be used to mea - 
sure the black hole mass quite precisely l|Rusin et al.ll200^ . 
The second central image would be a saddle rather than a 
maximum, so it would be interesting to consider whether 
a central saddle would have different microlensing proper- 
ties than the central maximum we have considered. We sus- 
pect that, at the radio wavelengths where this application 
would be pursued, microlensing fluctuations are still small 
compared with changes in the smooth lens model. Never- 
theless, it would be interesting to understand in detail how 
microlensing of central saddles compares with microlensing 
of central maxima. 
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